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ABSTRACT 



Context. The analysis and interpretation of the H 2 line emission from planetary nebulae have been done in the literature by assuming 
that the molecule survives only in regions where the hydrogen is neutral, as in photodissociation, neutral clumps, or shocked regions. 
However, there is strong observational and theoretical evidence that at least part of the H 2 emission is produced inside the ionized 
region of these objects. 

Aims. The aim of the present work is to calculate and analyze the infrared line emission of H 2 produced inside the ionized region of 
planetary nebulae using a one-dimensional photoionization code. 

Methods. The photoionization code Aangaba was improved in order to calculate the statistical population of the H 2 energy levels, as 
well as the intensity of the H 2 infrared emission lines in the physical conditions typical of planetary nebulae. A grid of models was 
obtained and the results then analyzed and compared with the observational data. 

Results. We show that the contribution of the ionized region to the H 2 line emission can be important, particularly in the case of 
nebulae with high-temperature central stars. This result explains why H 2 emission is more frequently observed in bipolar planetary 
nebulae (Gatley's rule), since this kind of object typically has hotter stars. Collisional excitation plays an important role in populating 
the rovibrational levels of the electronic ground state of H 2 molecules. Radiative mechanisms are also important, particularly for the 
upper vibrational levels. Formation pumping can have minor effects on the line intensities produced by de-excitation from very high 
rotational levels, especially in dense and dusty environments. We included the effect of the H 2 molecule on the thermal equilibrium of 
the gas, concluding that, in the ionized region, H? only contributes to the thermal equilibrium in the case of a very high temperature 
of the central star or a high dust-to-gas ratio, mainly through collisional de-excitation. 

Key words. Astrochemistry - Infrared: ISM - ISM: molecules - planetary nebulae: general - ISM: lines and bands 



1. Introduction 

Sinc e the first detectio n of H2 in a planetary nebula (PN) 
by ITreffers et al.l (Il976h . this molecule has been detected in 
many PNe both i n ultraviolet (UV) absorption and infrared 
(IR) emission (e.g. iBeckwith et al. 1978h IZuckerman & Gatleyl 



genitor dGussie & Pritchetlll988t iNatta & Hollenbachll 1998b. or 
neutral clumps inside the ionized region jBeckwith et al. 1978; 
i Gussie & Pritched l 19881: Ikeav et al.lll988t lTielenslll993h Ischildl 
ll995USpeck et alj|2002b . 



1988; Aspinetal. 1993; Kastneretal. 
McCandliss et al.N2007t ISterling et all 

2004). The H2 molecules can be excited by UV photons, colli- 



2005; Herald & Bianchi 



sions with the gas, or by formation on excited levels. The con- 
tribution of each mechanism to the population of the H2 lev- 
els depends on the physical conditions of the gas. Up to now, 
analyses of the excitation mechanism of the H2 mol ecules pro- 
ducing the observed IR li n es have been inconclusive ([Dinersteinl 



On the other hand, there is strong observational and the- 
19961 : iHora et alJ 1999 : oretical evidence that at least part of the H2 emission is pro- 
duced inside the ionized region of PNe. Some authors have 
noticed that the morphology of some PNe shown by images 
taken in the H2 1-0 S(l) infrared line is very similar to those 
taken in [Nil], [S II], and [OI] forbidden optical lines and 
in hydrogen re combination lines, which ar e produced in the 



clucing tne observed IK lines nave been inconclusive (Dinerstein 
199UlShupe etaD [1991 IHora et al] [19991 ISoeck et al J 12003: 
Rosado & Ariasll2003l:lLikkel et al.ll200 6: Matsu ura et al.l2007h 



In the literature, the H2 IR emission lines from PNe is usu- 
ally analyzed under the assumption that H2 molecules only exist 
in regions where the hydrogen is neutral, such as photodissoci- 
ation regions (PDRs; [Tielensll 1993t INatta & Hollenbachl [19981: 
IVicini et all [l 9991; iBernard-Salas & Tielensll2005l) . shocked re- 
gions between the expanding envelope and the wind of the pro- 



Webster et al. 


1988; Zuckerman & Gatley 


1988;' Balick et al. 


1991; Schildl 


1995; Allen etaljh 


9971 Guerrero et al. 20001: 


Lopez et al. 2000; Arias et al. 2001 


; Bohiaas 2001; Speck et al. 


2002, 120031). Excitation temperatures of approximately 1000 
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to 2000 K are i nferred from observations of Fb IR emission 
from some PNe (IHora et alJll999l iLikkel et alJl2006L and refer- 
ences therein). Such high excitation temperatures indicate that 
the molecule is excited by a strong UV radiation field, since 
collisions or shocks cannot explain the presence of lines from 
excited levels with vibrational num bers over three dHora et aTl 
119991: iBlack & van Dishoecld[l9"87h . Furthermore, in a previous 
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paper eman & Gruenwaldl2004 hereafter Paper I), we calcu- 
lated the density of H2 inside the ionized region of PNe, showing 
that H2 can survive in a partially ionized region with moder- 
ate temperature where neutral and ionized species coexist. Since 
such regions can be large when the temperature of the ionizing 
star is high, the ionized region can be a potential contributor to 
the total H2 line emission in some PNe. 

In the present paper we study the contribution of the ionized 
region to the observed H2 IR line emission of PNe. The H2 IR 
line intensities are calculated with a photoionization code. The 
effects of the temperature and luminosity of the central star, gas 
density, and dust-to-gas ratio on the line emission are studied. 
The models are described in Sect. [2] In particular, the formal- 
ism adopted for calculating the H2 energy level population and 
the intensities of the H2 emission lines, as well as the included 
H2 excitation mechanisms are described in detail. Results are 
discussed in Sect. [3] A summary of the conclusions and final 
comments are presented in Sect. [4] 

2. Models 

To calculate the intensity of H2 infrared lines emitted by the 
studied region, the molecular density and the population of each 
rovibrational energy level of the electronic ground state must be 
known. Both the molecular density and the level population de- 
pend on the location inside the nebula, since the radiation field, 
and thus the physical conditions of the gas, depends on the dis- 
tance to the ionizing source. The physical conditions of the nebu- 
lar gas were ob tained with the one-dimens ional photoionization 
code Aangaba (iGruenwald & Vieg as 1992). For a descripti on of 
the code and of the calculation of the H2 density see iPaper fl . 
For the present study, the calculation of the H2 level population 
and line intensities were introduced into the numerical code. The 
formalism adopted for these calculations is described below. 



2.1. The H 2 energy level population 

In Aangaba, the population of each H2 level can be calculated 
by assuming either local thermodynamical or statistical equi- 
librium. In the first case, the level population is set by the 
Boltzmann distribution and only depends on the gas tempera- 
ture. We assume statistical equilibrium, where the population of 
the energy levels is given by the balance between the rates of 
population and depopulation for each level. This is equivalent 
to stating that there is no net variation in the population of each 
level, i.e., 



dn w (y, J) 
dt 



0. 



(1) 



In the equation above, n w (v, J) is the density of H2 in an energy 
level with vibrational and rotational quantum numbers, respec- 
tively, v and J of the electronic state w. There is an equation as 
Eq. [T|for each level. This set of equations is linearly dependent 
and, to solve it, we replace one of the equations by the equation 
of conservation 



2n w (v,7) = n(H 2 ) 



(2) 



where n(H 2 ) is the density of molecular hydrogen. 

The expression on the left hand side of Eq. [T] is evaluated 
by the algebraic sum of the population and depopulation rates 
of the level (v, J) owing to the various mechanisms. To calculate 



the population in the rovibrational levels of the ground state, we 
included several excitation and de-excitation mechanisms (radia- 
tive and collisional), as well as the possibility that H2 is produced 
or destroyed at a given level by chemical reactions. Equation [1] 
can be rewritten as 



dn w (v, J) 



dt 



dn w (v, J) 



rad 



dt 



dn w (v, J) 



col 



dt 



= 



(3) 



chem 



where the terms on the left hand side of the equation are the 
terms that correspond to level population, respectively, by radia- 
tive transitions, collisional transitions, and formation or destruc- 
tion processes by chemical reactions. The processes included in 
the calculation are described in Sects. IXTTTI to 12. 131 To calcu- 
late the population of the upper electronic states, only radiative 
electronic transitions between the excited state and the ground 
electronic state are taken into account (see justification in Sect. 
12X21 

The set of equations can be solved for the density of each 
level if the rate coefficient for each process and the densities of 
the reactant species are known. The total H 2 density, as well as 
the densities of H°, H + , H~, Hj, and Hi, is given by the solution 
of the chemical equilibrium equations (see Paper I). In the fol- 
lowing we discuss the included mechanisms of populating and 
depopulating the energy levels of the H2 molecule. 

2.1.1. Radiative transitions 

The population of the H2 rovibrational levels of the electronic 
ground state by radiative mechanisms occurs through two main 
routes: electric quadrupole transitions between the rovibrational 
levels, involving IR photons, or electric dipole transitions to up- 
per electronic states with posterior decay to the ground state, 
involving UV photons. 

The transitions between two energy levels with different ro- 
tational and/or vibrational quantum numbers within the same 
electronic state are called rovibrational transitions. For the elec- 
tronic ground state of H2, X'E+ (hereafter only referred to as X), 
they occur through electric quadrupole transitions, since transi- 
tions by electric and magnetic dipole are forbidden. Quadrupole 
transitions have much lower probabilities than the electric dipole 
ones. For the X state, we included 289 rovibrational levels dis- 
tributed over 15 vibrational levels. T he H 2 rovibrationa l transi - 
tion probabilities were calculated by IWolniewicz et ail (Q998), 
and the energies of the levels were kindly provided by E. Roueff 
(2005, private communication). The selection rules for these 
transitions imply that the change in the rotation quantum num- 
ber must be AJ - 0, +2. The lines produced by rovibrational 
transitions of H2 are in the 0.28 /im to 6.2 mm range of the elec- 
tromagnetic spectrum. 

Transitions where the electronic state of the molecule 
changes, with or without change in the rotational or vibrational 
quantum numbers, are called electronic transitions. Electronic 
transitions are allowed by electric dipole and are thus very 
likely to occur. They are efficient in exciting and destroying H2 
molecules and must then be included in the calculation of the sta- 
tistical population of level X. We included electronic transitions 
between the ground state and the excited electronic levels B'E+ 

and C'lTu (Lyman and Werner bands, respectively). The levels 
B'E+ and C l Yl u are referred to hereafter just as B and C, respec- 
tively. The A doubling splits the energy levels of C state in two, 
denoted here by C and C + . For the B electronic level, 408 rovi- 
brational levels, distributed over 38 vibrational levels, are taken 
into account; for C and C + , 143 and 142 rovibrational levels, 
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respectively, are included, distributed over 14 vibrational levels. 
Excitation to B and C levels can be followed by decay back to 
the ground state on any rovibrational level (this route is known 
as UV pumping) or to the vibrational continuum of X leading to 
the molecule dissociation. This last route is called photodissoci- 
ation in two steps or Solomon process, and it is described in Sect. 
12.1.31 The Einstein coefficients and the wave number of the elec- 
tronic transitions have been calculated by Ab g rail et al.l 1 1994) 
and the photodissociation fractions by Abgrall et al. (2000). The 
selection rules for electronic transitions requires that the change 
in the rotation quantum number must be AJ = +1 for Lyman 
band transitions, while in the case of Werner band transitions 
AJ = for X-C transitions and AJ = +1 for X-C + transi- 
tions. The electronic excitation of the Lyman and Werner bands 
requires photons with energy between 6.7 and 15 eV. 

2.1.2. Collisional transitions 

Collisions may also change the H2 energy level. In photoion- 
ized regions the average energy of the particles is less than 1.7 
eV, and in regions where the H2 density is significant the aver- 
age energy is even less. Since high energies are needed for elec- 
tronic transitions (> 6.7 eV), collisional electronic excitation is 
unlikely to occur. On the other hand, collisional rovibrational 
transitions may be very important in astrophysical environments 
and are so taken into account. We included collisions of H2 with 
the main components of the gas, i.e., H, H + , He, H2, and elec- 
trons. Collisions may change the total nuclear spin by inducing a 
temporary non-zero dipole moment in H2 (the so called reactive 
collisions, in opposition to the non-reactive collisions, in which 
there is no change in the molecular nuclear spin). 

Rate coefficients for the de-e xcitation of H? by coll isions 
with H atoms were calculated by iMartin & Mandvl ( 1 19951) . us- 
ing the quasiclassical trajectory method, for all H2 rovibrational 
levels of the ground state, for both reactive and non-reactive col- 
lisions. l^jiourlotetaD ([l999|) give a compilation of references 
for the rate coefficient of non-reactive collisional de-excitation of 
H2 by H, He, and H2 calculated by fully quantum methods, but 
only for H2 vibrational levels up to v = 3. This method of calcu- 
lation is more accurate than the quasiclassical trajectory method 
(D. Flower 2005, private communication) and was then preferred 
when available. 

Rate coefficients f or de-ex c itation of H2 by collisions with 
H + were calculated bv lGerlic h ( 1990) for pure rotational transi- 
tions of v = 0, with J upper up to 9 and A J up to 9. The remaining 
pure rotational transitions of v = (J upper from 10 to 29) were 
estimated by extrapolating the existing coefficients. To approx- 
imately reproduce the behavior of the available coefficients, we 
assume the following extrapolations: 

- For Ji ower < 7, we assume the same coefficient of the transi- 
tion from the h ighest J upper of the same parity calculated by 
lGerlichl (fi990). going to the same lower J level of the tran- 
sition in question. For example, the coefficient for the transi- 
tion from J U pp er — 1 1 to Jiower — 2 is equal to the coefficient 
for the transition from J upper = 9 to Ji ower = 2. 

- For transitions with J[ ower > 8, we assume t he same coeffi- 
cient of the highest transition calculated by iGerlichl d 19901) 
involving upper and lower J levels of the same parity. For 
example, the coefficient for the transition from J upper = 14 
to Jhwer — 9 is equal to the coefficient for the transition from 

Jupper — 8 to Jiower ~ 7. 

For pure rotational transitions of levels v >1, we assume 
the same coefficients as the corresponding transition of v = 0. 



Coefficient rates for H2 vib rational de-excitati on by electrons 
for Av up to 3 are given by Draine et al. ( 1983). Since they do 
not provide the /-resolved rate coefficients, we assume that the 
transition rate for each rovibrational transition is proportional to 
the Einstein coefficient of the corresponding radiative transition. 
Only non-reactive transitions are included {AJ = 0, +2). 

In the present work we assume that the coefficients for the 
collisional excitation of H2 by electrons are the same as those 
for collisions with H + , because of the lack of rate coefficients 
for this mechanism. 

2.1 .3. Formation and destruction of H 2 at a given energy 
level 

In iPaper J the main processes of formation and destruction of 
H2 due to chemical reactions in the ionized region of PNe were 
determined. Chemical reactions form or destroy H2 in a given 
level. Since the H2 molecule can be significantly produced or 
destroyed by these mechanisms, they must be included in the 
set of equations of statistical equilibrium. Unfortunately, level- 
resolved rate coefficients for most of the reactions are still poorly 
known and some assumptions must be made. The processes and 
our assumptions are discussed in the following paragraphs. 

Photodissociation - There are two main routes for H2 pho- 
todissociation: the direct and the two-step processes. In the di- 
rect process, the H2 molecule is excited from the ground state 
to the vibrational continuum of an upper electronic state by an 
UV photon. The cro s s se ction for this process was calculated by 
Allison & Dalgarnoi d!969t) for each of the 15 bound vibrational 
levels of the ground state of H2. However, since they provide no 
information on the cross section for each rotational level, we as- 
sume that the coefficient rate is the same for all rotational levels 
of the same vibrational level. 

As said above, the H2 photodissociation in two steps be- 
gins with the excitation of the molecule from the ground elec- 
tronic state to a bound rovibrational level of an upper elec- 
tronic state. The subsequent decay to the vibrational contin- 
uum of the ground electronic level leads to the d issoc iation 
of the molecule. Ac cording to IStecher & Williams ( 19671) and 
lAbgrall et al.l {1997), this happens in approximately 1 1% of the 
excitations to the B and C levels. This process is particularly im- 
portant because molecules can be dissociated by photons with 
energies smaller than the H ionization potential. As a result, the 
rate of this process is not significantly affected by the H column 
density, but depends mostly on the geometrical dilution of the 
UV radiation and on the H2 column density (i.e., depends on the 
self-shielding). The Einstein coefficients for the electronic ex- 
citations and the fraction of the molecules t hat dissociates (dis- 
socia tion fraction) have been obtained by Abg rall et al.l d 19941 
2000), respectively. The effect of the H2 self-shielding on the 
photo dissociation in two step s was included, using the formal- 
ism of Bl ack & van Dishoeckl(ll987l) . 

Photoionization - The threshold energy for the photoioniza- 
tion of H2 is £^(0,0) = 15.4 eV for the ground rovibrational 
level of X and is smaller for upper rovibrational levels. It is as- 
sumed here that all molecular ions are formed at the ground 
level. For H2 in vibrational levels v > 4, the threshold energy 
of photoionization is lower than the photoionization potential 
of H (13.6 eV), but since the population of such levels are of- 
ten very small, particularly in the recombination zone where the 
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H 2 density is significant, this is not an important depopulation 
process. An accurate cross section for the photo ionization from 
the vi brational ground level of H2 is given by lYan et al] (119981 
I2001IP1 . We assume the same shape for the cross section for the 
upper vibrational levels, except that it is displaced in energy by 
the difference in the photoionization threshold energies of the 
levels. Since there are no /-resolved cross sections, it is further 
assumed that the photoionization rate of a given vibrational level 
is distributed uniformly among its rotational levels. 



assume that the H2 mole cule is dissociated fo llowing the level 
distribution suggested bv lLenp&Shul(ll983l) . 



Formation of H 2 on grain surfaces - In lPaper fl, we described 
the model adopted for the grain surface reaction and the expres- 
sion used for its rate coefficient. The coefficients for H2 forma- 
tion by this rea ction for each energy level have been disc ussed 
in the literature dTakahashi & Ueharal200lllTine et alj2003l and 
references therein), but they are not well established yet. We as- 
sumed the sam e level distribution of H2 produced by this process 
as adopted bv lBlack & Dalgarnoldl976l) . 



Associative detachment reaction 

ment, 

H + H^H 2 + hv, 



The associative detach- 



(4) 



is a major route for the f ormation of H2 in the ionized region of 
PNe (Blackl fl978l iPaperll) . The rate coefficients for the forma- 
tion of H2 at each b ound level of t he gro und state by this reaction 
were calculated bv lLaunav et alj d!99 lh . 

Charge exchange reactions - The H 2 molecule can be pro- 
duced and destroyed by charge exchange reactions 



H + + H 2 



H + + H 2 



H + H+. 



(5) 



(6) 



The rate coefficie nts for these p r ocesse s are given by 
iKarpas et alJ (Il979h and iGalli & Pailal (1 19981) . respectively. For 
the first process we assume that the level distribution of the 
res ulting molecules has the s ame distribution as the one used 
by iBlack & Dalgarnd ([1976) for molecules formed on grain 
surfa ces, but with the app ropriate energy of formation (1.83 
eV, iHollen bach & McKeelTT979h . For the destruction route we 
adopted the formula for the coefficients given bv lGam&Paijp 
(1998). Since the level-resolved coefficient rates are unknowrfl, 
we assume that H2 molecules can be destroyed by charge ex- 
change reaction in any rovibrational level equally. 



Collisional dissociation - For the collisiona l dissociation of 
H? by H atoms we used the formula given by IShapiro & Kangl 
(198 7]) for the rate coef ficient, while the coefficient is that of 
Don ahue & Shulll (1 1 99 lh for the dissociation by electrons. We 



1 There is a small correction to their fitting formula for the first en- 
ergy interval: the coefficient of the second term should be 197.25 (H. 
Sadeghpour 2010, private communication). 

After the conclusion of this work we found that lSavinetaLl d2004h 
had published coefficient rates for the H 2 destruction in each vibrational 
level by this process. However, as we show in the following, this reac- 
tion (and its reverse) is not an important mechanism of H 2 excitation, 
so our conclusions will not be affected by using these new coefficients 



k(v,J) = k T0JAh - 



E(y,J) 



E(v,J) 



(7) 



where &total is the total coefficient rate of collisional dissocia- 
tion of H2, and E(v, J) is the energy of the level (v, J). 

2.2. H 2 and the thermal equilibrium 

Thermal equilibrium is assumed for the calculation of the gas 
temperature, that is, the total input of energy in the gas per unit 
time and volume is balanced by the total loss of energy per unit 
volume. Besides the processes of gas heating and cooling asso- 
ciated with the atomic species and with the g rains (for details 
of ho w dust is included in the calculations, see lGruenwald et ail 
2011), in the present work, we also included the H2 molecule 
contribution to the energy balance of the gas. The mechanisms of 
gas heating due to H2 are collisional de-excitation, photoioniza- 
tion, photodissociation, formation on grain surfaces, associative 
detachment, and charge exchange reaction; the mechanisms of 
energy loss are collisional excitation, charge exchange reaction, 
and collisional dissociation. 

2.3. H 2 infrared emission line intensities 

The emissivity of a line produced by de-excitation from a given 
level is the product of the population of that level by the Einstein 
coefficient of the transition and by the energy of the emitted pho- 
ton. The emissivities of the H 2 IR lines are calculated for each 
point of the nebula. In our models, we assume that the nebula 
is radiation-bounded. The intensities are calculated by integrat- 
ing the emissivity over the nebular volume. It is also possible 
to obtain solutions for matter-bounded PNe by integrating the 
emissivity up to a given external radius. The gas is assumed to 
be optically thin for the H2 lines. 

2.4. Input parameters 

A grid of theoretical models was obtained for incident radi- 
ation spectrum, gas density, and dust-to-gas typical of PNe. 
The central star is assumed to emit as a blackbody. The gas 
density and the elemental abundances are assumed constant 
along the nebula. We use the mean values for PNe obtained by 
Kingsburgh & Barlowl d 19941) for the abundances of He, C, N, 
O, Ne, S, an d Ar. For Mg, Si, CI, and Fe, w hose abundances are 
not given by Kingsburgh & Barlow ( 1994), the values are those 
adopted by Stasihs ka & Tvlenda (119861) . in order to roughly cor- 
rect for grain depletion. Since the effect of these elements is not 
very important for the gas cooling, their exact proportion in the 
form of grains is not important for our results. Amorphous car- 
bon grains with Q.l/im radius are assumed. 

We define a standard PN with a given set of input parameters 
to explore their effect on the level population and on the inten- 
sity of the lines emitted by the H2 molecule. For this, we vary 
one of the parameters within its typical range, while keeping the 
other parameters fixed. Unless otherwise noted, the parameters 
of the models are those of the standard PN. The input param- 
eters for the standard PN are the following: = 150000 K, 

= 3OOOL0, «h = 10 3 cm -3 , and the assumed dust-to-gas 
ratio is M^/M g = 10~ 3 . The temperature chosen for the central 
star of the standard PN is above the average for observed PNe (~ 
80 000 K) in order to highlight the effects of H 2 . 
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Fig. 1. The standard PN model, (a)-(c) Ionizing continuum spec- 
tra in three different positions of the nebula, as indicated in the 
cartoon representing the ionization structure of the standard PN 
in panel (d) (the PDR is not to scale), (e) Radial density profile 
of H°, H + , H~, H2, , and Hi in the ionized region. Densities 
are relative to the total density of H nuclei, (f) Radial gas temper- 
ature profile in the ionized region. The blue band in both panels 
indicates the TZ. A color version of this figure is available on- 
line. 

3. Results 

3.1. The transition zone 

In our discussion we formally define the outer boundary of the 
ionized region as the location where one of the following con- 
ditions is reacheofl a) the fractional abundance of ionized hy- 
drogen reaches 0.01% or b) the gas temperature is l ess than 
100 K . The ionization degree value is chosen following lTielensI 
(2005). Since for some PNe models this ionization degree is only 
reached in very low gas temperatures, we include the additional 
temperature stop criterion. In most of our models, the latter con- 
dition is reached first. We also define the transition zone (TZ) 
as the region between the location where the fraction of ionized 
hydrogen is 95% and the outer boundary of the ionized region. 
A cartoon representing the ionization structure of a PN is shown 
in Figure [TJl. 

Figure[T]shows results for the standard PN. The ionizing con- 
tinuum spectrum is given for three different positions in panels 

3 We note that the same criteria were assumed in iPaper j models, 
although, due to an oversight, only the ionization degree was mentioned 
in that paper. 



(a) to (c). Panel (c) shows the incident continuum spectrum at 
the inner edge of the NP. The spectrum is the geometrically di- 
luted stellar blackbody. Panel (b) shows the continuum spectra at 
the inner edge of the TZ. The spectrum shows absorption mainly 
due to hydrogen and helium, as well as diffuse UV emission fea- 
tures, but there is a significant flux of UV and soft X-ray pho- 
tons. These high-energy photons emitted by the central star are 
responsible for the formation of the TZ and are very important 
for the physics and chemistry in this zone. Panel (a) shows the 
continuum spectrum emerging from the ionized region. The flux 
of photons with energy above 13.6 eV is significantly absorbed, 
indicating the inner boundary of the PDR (see Tielens (2005). 

The radial distribution of H2 density (relative to the total H 
nuclei density) for the standard PN model is shown in Fig. QJ. 
In the figure, R is the distance to the central star. Relative densi- 
ties of H°, H + , H , Hi, and H^ are also included in this figure. 
The gas temperature profile is shown in Fig.QJ. The TZ for this 
model is the region indicated by the blue band in Figs. QJ and 
f. From these figures, it can be noticed that the relative density 
of H2 molecules reaches a m aximum in this warm and partially 
ionized region, as shown in PaperJ. For typical PNe models, 
the maximum relative H2 density in the ionized region is around 
10~ 5 , with the exception of models with hot central stars (T+ > 
250 000 K) where the density can be orders of magnitude higher. 
After a slight decrease, the rising of the H2 density to wards th e 
more neutral zone depends on the grain density (see IPaper J) , 
since the formation of H2 on grain surfaces is the most impor- 
tant processes in this outer zone of the ionized region. 

High-energy photons (hv > 100 eV) emitted by the central 
star are responsible for the formation of the TZ. They can pen- 
etrate deep into the nebulae, producing this extended region of 
partially ionized and warm gas. Such physical conditions favor 
the formation and survival of H2, as showed in lPaper J . The for- 
mation rates of H2 molecules due to associative detachment of H 
and H and to charge exchange reactions between H and H£ are 
enhanced in a mild temperature gas where neutral and ionized 
species coexist. Furthermore, the destruction is not very effec- 
tive because of self-shielding and shielding provided by atoms 
and grains present internally. Photons in the 13.6 to 100 eV range 
are absorbed by H and He atoms in the inner regions and do not 
contribute to the destruction of H2 molecules in this region. The 
amount of H2 inside the ionized region of PNe is then very sen- 
sitive to the thickness of the TZ. 

The available number of high-energy photons increases 
strongly with T*. As a consequence, the amount of H2 in the 
ionized region of PNe also increases, as can be seen in Fig. [2] 
where the H2 to total H mass ratio (Rm) is shown as a function 
of T+. Curves for different values of L*, tin, and Md/M g are plot- 
ted, respectively, in Figs. 2a, 2b, and 2c. 

The size of the TZ does not depend significantly on the lu- 
minosity of the central star. However, a luminous central star 
produces more photons and is more effective in photoionizing 
and photodissociating molecular hydrogen. Consequently Rm 
decreases for models with higher L* (Fig.|2^). 

The value of Rm also increases significantly with decreasing 
total H nuclei density («h; see Fig. \2)p). This effect is stronger 
in models with hotter central stars. The stronger absorption of 
the radiation in the gas in models with high density produces 
thinner TZs. It is important to notice that we are also increasing 
the dust density when we increase the gas density in a model, 
since the dust-to-gas ratio is fixed in any given model. The effect 
of changing the dust-to-gas ratio (Md/M g ) on the size of the TZ 
is complex since dust modifies the temperature of the gas, ab- 
sorbs the UV radiation, and acts as a catalyst for the molecular 
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Fig. 2. H2 to total H mass ratio {Rm) as a function of the stellar 
temperature. Curves for different (a) L*, (b) and (c) Md/M g 
are shown. The standard PN values are adopted for the parame- 
ters not mentioned. 



formation. As a result, the size of the TZ increases with Md/M g 
for low-temperature central stars, while the opposite effect oc- 
curs with hot stars. However, Rm increases when Md/M g is in- 
creased in the models, as can be seen in Fig. |2j;. This is mostly 
due to the increase in the rate of the reaction of H2 formation 
on grain surfaces. When significant density of grains is available 



(i.e., Md/Mg > 10 3 ), the reaction of H2 formation on grain sur- 
faces is an important process of H2 production in the outer part 
oftheTZ. 

Self-shielding is included in the present work (Sect. 2), re- 
su lting in an enhancement in Rm in comparison with the models 
of Paper I, particularly for models with high central star temper- 
ature, low central star luminosity, low gas density, and/or high 
Md/M g (compare Fig. |2] with their Fig. 4). 

3.2. The effect of H 2 on the thermal equilibrium 

Atomic species dominate the heating and cooling of the gas in 
the ionized region of PNe. If grains are present, they contribute 
significantly to the gas temperature (heating the gas) in a narrow 
zone near the inner border of the nebula and in the TZ. The in- 
crease in the temperature may be several thousand Kelvin in the 
first case and of a few hundred in the second. 

The effect of H2 on the thermal equilibrium is insignificant in 
most models, except in the cases of very high central star tem- 
perature and/or dust density. In these models, the net effect of 
the molecule is the cooling of the gas in the outer zone of the 
TZ, primarily through collisional excitation. For PNe with = 
350 000 K, for example, the gas temperature is decreased by up 
to 5000 K, while for the standard PN the temperature structure 
is not altered by the presence of H2. 

3.3. The H 2 1-0 S(1) line emission 

The 1-0 S(l) line is produced by the radiative transition from 
the level (v, /) = (1,3) to the level (0,1) of the ground electronic 
state and is one of the most intense lines of H2. Its wavelength 
is 2.122 /im, which is in the K atmospheric window, so it can be 
detected by both ground- and space-based t elescopes. This lin e 
has been detected in many PNe (for example lHora et ai] [T999). 

The fraction of H2 molecules at the level (1,3) as a function 
of the distance to central star is shown in Fig. [3] for the standard 
PN. This fraction varies by a factor of 5, while the H2 density 
varies more than 10 orders of magnitude (Fig. [TJ. For all our 
models, the density of H2 molecules in the level (1,3) and, as a 
consequence, the emissivity of the 1-0 S(l) line follow a radial 
profile similar to «(H2). Figure|4]shows the emissivity of the line 
H2 1-0 S(l) in the ionized region as a function of the distance 
to the central star for the standard PN model. The similarity of 
the radial profiles of the H2 density and the 1-0 S(l) emissivity 
can be noticed by comparing Figs. Q] and [4] The emissivity of 
this line is more important in the TZ, making this zone the most 
important contributor to the total 1-0 S(l) line intensity of the 
ionized region. 

According to our calculations, the level (1,3) is populated by 
both collisional and radiative mechanisms. The contribution of 
the chemical processes is negligible. For the standard PN, for ex- 
ample, collisions dominate the entire ionized region, while UV 
pumping is significant in the TZ. The relative importance of col- 
lisional over radiative processes for the level (1,3) population 
increases with increasing and tin, and also for decreasing L*, 
although this last effect is not very significant. Collisions domi- 
nate the excitation for models with > 200 000 K. Formation 
of H2 on grain surfaces may have some effect on the H2 energy 
level distribution, particularly for the case of higher rovibrational 
levels and high grain-to-gas ratio, but for the level (1,3) this ef- 
fect is not important. Other processes of H2 formation and de- 
struction of H2 do not contribute significantly to the energy level 
distribution of the molecule in the ionized region. 
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Fig. 3. Relative population of some H2 rovibrational energy lev- 
els of the ground electronic state as a function of the distance 
to the central star for the standard PN model. Levels are indi- 
cated by their vibrational and rotation numbers (v, J). The gray 
band indicates the TZ. A color version of this figure is available 
online. 
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Bry, [S II], [Nil], and [OI] as a function of the distance to the 
central star for the standard PN model. The gray band indicates 
the TZ. A color version of this figure is available online. 



In addition to the 1-0 S(l) line emissivity, Fig.|4]also shows 
the radial profiles of the emissivity of some atomic lines as cal- 
culated by our code. The lines are Bry in 2. 166 /an, [N II] /i6583, 
[O I] A6300, and [S II] AA61 1 6, 673 1 . These lines are mainly pro- 
duced in the region where H is ionized. According to our models, 
the peak in the emissivity of the 1-0 S(l) line in the ionized re- 
gion is very close to the region where the atomic lines [Nil], 
[OI], and [S II] are mainly produced. 



As mentioned in the Introduction, similarities between the 
morphology of images in the 1-0 S(l) line and those taken in the 
[Nil], [S II] and [OI] forbidden optical lines and hydrogen re- 
combination lines (H/3, Ha, and Bry) have been noticed in some 
PNe. Since the H2 1-0 S(l) line can be produced in a more ex- 
tended region than those lines, the similarity in the images is 
evidence that a neutral envelope is absent in those objects. In 
this case, all the H2 emission is produced in the ionized region, 
and the molecule coexists with N + , O , and H + in the TZ. 

A correlation between t he intensities of th e 1-0 S(l) and [O I] 
lines in PNe was found bv lReav et al.ldl988l) . Although they as- 
sociate this correlation with the existence of clumps inside the 
ionized region, it can also be naturally explained by both lines 
being produced in the same region, since the production of the 
[OI] line is very effective in the TZ, as shown in Fig. [4] 

Since the 1-0 S(l) line is mostly produced in the TZ, the 
flux of this line is very sensitive to T+, as can be seen in the 
three top panels of Fig. [5] where the flux of the 1-0 S(l) line is 
shown as a function of T+. Curves for different (left panels), 
n H (middle), and M d /M g (right) are shown. The 1-0 S(l) line flux 
does not depend much on L + or Md/M g , but this flux decreases 
strongly if the model has a gas density higher than the stan- 
dard PN value of 1000 cm -3 . The strong increase in the 1-0 S(l) 
flux with T+ we obtai n with our mod els agrees with the results 
and conclusions from Phillips (2006). After analyzing observa- 
tional data from the literature, Phillips (2006) find a correlation 
between the detection of H2 in PNe and their Zanstra tempera- 
tures, which he attributes to an effect of soft X-rays. The strong 
increase in the 1-0 S(l) flux with can also explain why the 
detection of H2 is more common in bipolar PNe (Gatley's rule; 
IZuckerman & Gatlev|[l9 88; Kast ner et al.lll996|). sin ce, accord- 
ing to lCorradi & Schwarzld 19951) andlPhillipsI d2003l). th ese ob- 
jects typically have higher T+. [StanghellinietaLl d2002) do not 
find a similar correlation in their analysis. This, however, may be 
due to their selection criterion that excludes PNe with high un- 
certainties in Zanstra temperature and, therefore, excludes high 
Zanstra temperature objects. 

The ratio of the total luminosity in H2 IR lines to the total 
luminosity of the central star is shown in the panels of the second 
row in Fig. [5] The sensitivity of this ratio with T* is evident. PNe 
with smaller L+ or «h also convert the stellar continuum in H2 
emission more efficiently. The change in the dust to gas ratio 
have a small effect in this ratio. Note that the total luminosity 
emitted by the H2 lines may reach up to a few percent of the 
luminosity of the central star for models with a very hot and 
low-luminosity central star. This ratio is much lower for more 
typical PNe. 

The 1-0 S(l) line intensity is around 1% of the total H 2 IR 
emission in most models, as shown in the plots of the bottom 
row in Fig. [5] The ratio of the 1-0 S(l) line intensity to the total 
H2 IR line emission does not depend much on L* or Md/M g , but 
it increases for models with higher and decreases for models 
with «h > 10 3 cm 4 . 

The ratio 1-0 S(l)/Bry in the ionized region as a function 
of the PNe model parameters is shown in Fig. [6] Models 
with different L+, tin, and M^/M g are shown in panels a, b, 
and c, respectively. Ratios obtained from observations are 
also shown . The observe d value s are from the f ol lowing 



references: Beckwith et a il d 19781). Ilsaacmanl (fl98l. IStorevi 
(1984), Webst er et al.l d 19881). lAspin et al.l d 19931). lAllen et al ' 
A 1997k Mora et al l (1 19991). [Davis et all (l2003h. iGuerrero et a! 
(hood). iRudvetalJ <l200lk . iKellv & Hrivnakl d2005l) and 
Likk el et al.l (120061) . In this fi gure we include ratios obtained 
from narrowband images and from long-slit spectra. In the latter 
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Fig. 5. Flux of the 1-0 S(l) line in arbitrary units (top row panels), ratio of the total H2 IR lines luminosity to the central star 
luminosity (middle), and 1-0 S(l) luminosity to total H2 IR lines luminosity ratio (bottom) as a function of r*. Curves for different 
L+, rin, and Md/M g are shown, respectively, in the left, middle, and right panels in each row. The line styles are the same as in Fig. 
[2] In the bottom left panel, the curves for L* = 100 and 3000 Lq coincide. The standard PN values are adopted for the parameters 
not mentioned. 



case, the slit is often positioned here the H2 emission is more 
intense, so the ratio 1-0 S(l) to Bry may be biased to higher 
values. In Fig. [6] represents the temperature of the stellar 
blackbody for models and the Zanstra or the energy-balance 
temperature for observations. For the Zanstra temperature, 
we used Tz(Hell) when available and Tz(HI) otherwise. In 
the few cases we did not find the Zanstra temperature, we 
adopted the value calculated throug h the energy-balance method 
dPreite-Martine z & Pottasch 1983). When more than one value 
(obtained by the same method) is found, an average is assumed. 
The temperatures we r e taken from th e following r e ferences: 
iPottasch et all (fl978h . iMartinl jfl98lb. iKalerl ( fl98l IPottascT 



(119841). iReavet all (11984. IShaw & Kalerl (I1985L 1198 
de Freitas Pacheco et al l (1 1 9861). iGathier & Pottasch! 119881 
1989 ). iGleizes et : al.l dl989l). IKaler & Jacobvl (11989. fl99ll) 



1989), Gleizesetal. (1989), Kaler & Jacobv (1989 
Jacobv & Kalerl (Il989h. iPreite-Martinez et all dl989 



Walton et al.1 d 19891). IKaler et al.1 (1 19901). iMendez et al.l ( 
Stasinskaet al J (119971) . iBohigasI d2001l) , and~~ 
d2009l) . 



1991 



), 

1992b 
Szvsz ka et alj 



Planetary nebulae can be radiation- or matter-bounded, de- 
pending on the amount of matter in the nebula and on the ion- 
izing spectrum of the central star. In the first case, the ionizing 
spectrum cannot ionize the whole nebula. As a result, there is 
a neutral outer region. The emission of Bry is produced inside 



the H ionized region, while the emission of 1-0 S(l) line is pro- 
duced in both the TZ and the neutral region. In our models, we 
assume that the nebula is radiation-bounded, therefore the line 
intensities are obtained by integrating the emissivity along the 
nebula, from the inner to the outer radius of the ionized region 
of the PN. Since in this work only the contribution of the ionized 
region is calculated, in the case of a radiation-bounded PN, only 
a lower limit is obtained for the H2 line intensities and for the 
line ratio H2 to Bry. For a matter-bounded PN, the whole nebula 
is ionized, and our ratios are upper limits. 

Figure|6]shows that the ionized region can be responsible for 
part or even the whole 1-0 S(l) emission of PNe. The ionized 
region can account for the whole emission for some PNe with 
high T+. In these cases, the hot star may be ionizing the whole 
(i.e., the PN is matter-bounded) or most of the nebula. 

Most PNe with detected H2 emission, however, may 
have a molecular envelope, which could be responsible for 
the differences between the calculated and observed ratios. 
Contamination by other nearby H2 line s is negligible acc ording 
to our models. The spectral resolution in lHora et al.l (Q999), from 
which most of the observation data were obtained, allows resolv- 
ing the 1-0 S(l) line from a neighboring helium line, which may 
otherwise contaminate the molecular line. 
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Fig. 6. Line intensity ratio of H2 1-0 S(l) to Bry for the ionized 
region as a function of T*. Curves for different (a) L*, (b) riu, 
and (c) Md/Mg are shown. The line styles are the same as in 
Fig. 12 Observed ratios are represented by dots. Boxes with up 
arrows are lower limits and with down arrows are upper limits. 
References for the observations are given in the text. The stan- 
dard PN values are adopted for the parameters not mentioned. 



3.4. Other important H 2 lines 

The H2 molecule emits a few thousand rovibrational lines in the 
IR. According to our models, the most intense lines are emit- 
ted in the range 1 to 29 fim, most of them with wavelengths 
from 1 to 5 /mi (near IR). Above 5 fim the most intense lines are 
in general pure rotational lines, particularly S(0) to S(8) of the 
0-0 band. The majority of the intense lines are produced by 01- 
tho levels, especially those from higher vibrational levels, since 
the ortho-para ratio of H2 is normally higher than unity. Table 
Q] lists those lines that are more intense than 1 % of Bry in at 
least 10% of the models we ran, covering a variety of PN pa- 
rameters. Assuming this criterion, the most intense lines in the 
ionized region are lines of the molecular bands 0-0, 1-0, 1-1, 2-0, 
2-1, 3-1,4-1, 4-2, and 5-3. Several of these intense lines have al- 
ready been detected in PNe (they are indicated by an asterisk in 
Table [TJ. Se veral authors published observations of H2 lines in 
the K band (iRamsav et alJI 19931 iHora et all [19991: IVicini et all 
1999, for example). IHora et all ([1999) report the detection of 
more than 50 H2 lines, with upper vibrational levels with quan- 
tum number up to 11, in the J, H, and K bands, which are 
rich in H2 lines. They are atmospheric windows, which allows 
both ground and space-based observations. Molecular hydrogen 
lines were also detected in other s pectral ranges in PNe. For in- 
stance. lMatsuura & Ziiistral (l2~005) detected the lines S(l) to S(6) 
of the molecular band 0-0 in observa tions of NGC 6302 with 
the In frared Space Observatory (ISO); Bernard-Sala s & Tielensl 
(120051) published the detection of the lines 0-0 S(l) to S(7), 1-0 
Q(l), 1-0 Q(3), 1-0 Q(5), 1-0 0(3), 1-0 0(4), 1-0 0(5), 1-0 0(6), 
and 2-1 0(3), also in observations made with ISO; the line s S(2) 
to S(7) of the 0-0 band were detected bv lHora et ail (120061) in the 
Helix nebula (NGC 7293) with the Spitzer Space Telescope. 

Ratios of some H2 lines to Bry are shown in Fig. [7] The plots 
are similar to those shown in Fig. [6] As previously mentioned, 
in the ionized region, the production of the 1-0 S(l) line is more 
efficient in the TZ. Similarly, all H2 lines are produced mostly 
in this zone. As a consequence, the intensity of the H2 lines de- 
pends on the PN parameter in a similar same way as the 1-0 S(l) 
line, as can be seen in Fig. [7] 

Ratios obtained from observations of PNe are also included 
in Fig. [7] References for the central star tempe rature are the same 
as in Fig. [6] Th e line ratio s are ob tained from lHora et al.l (fl999) 
and [Sterling & Dinersteinl ([2008). Our result show that the ion- 
ized region can account for part or even all the H2 emission of 
these lines, particularly for PNe with hot central stars. 



3.5. Excitation mechanisms of the H 2 levels 

In all the analyzed models, collisional transitions are a major 
mechanism for population and depopulation of the H2 energy 
levels, in the whole ionized region. This mechanism dominates 
the excitation of the lower vibrational levels. Collisions with 
ionized species are important in the inner (hotter and most ion- 
ized) zone of the models. The net effect of this mechanism is to 
populate levels with J < 5 by de-excitation of the upper levels 
(especially with J > 7). Collisions with neutral species, on the 
other hand, are more important in the TZ, where the temperature 
is moderate and neutral species are abundant. The net effect of 
this mechanism is to de-excitate the levels with lower rotational 
number, populating the upper rotational levels, especially those 
with lower vibrational numbers. The importance of collisions in- 
crease with T+, since the TZ is more extended in this case. 

The rate of H2 rovibrational radiative excitation is negligible. 
On the other hand, rovibrational de-excitation following elec- 
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Table 1. Important H2 lines of the ionized region of PNe 



Line^ A(p.m) Line AQim) Line A(p.m) 



" Lines with (*) are already or possibly detected in PNe. By possibly 
detection we mean that a blend with or contamination by other line may 
prevent confirmation of the detection so far. 



tronic (UV pumping) or rovibrational collisional excitation is 
very important. The main effect of this mechanism is the ex- 
citation of H2 from the lower rotational levels of v = to the 
higher vibrational levels (v > 3) with J < 7, For levels with 
v > 9 and J < 7 UV pumping is the dominant mechanism. The 
contribution of the C and C + states for the UV pumping rate is 
similar to the contribution of the B state. UV pumping is partic- 
ularly important for H2 excitation in the case of PNe with lower 
«h and/or . 



Formation pumping has a secondary effect on the H2 popu- 
lation in the ionized region. Such an effect is only noticed for 
the upper levels (v > 5 and J > 6) in the outer fraction of PNe, 
especially for denser nebula or colder central stars. In this case, 
grain surface reaction and associative detachment are the main 
processes. Grain surface reaction may have a stronger effect for 
high dust density models. 

The H2 destruction processes do not affect the H2 population 
in the ionized region. Photoionization has a small effect (less 
than 10%) in populating of the lower vibrational levels in the 
inner zone of the nebula. However, since the H2 density is very 
low in this region, it does not affect the intensities of the H2 lines 
of PNe. 

4. Summary 

In the present work we have studied the H2 infrared emission 
from the ionized region of PNe. For this, we used the one- 
dimensional photoionization code Aangaba, in which we in- 
cluded the physics and chemistry of the H2 molecule. This pow- 
erful tool can now be used to study the H2 density, level pop- 
ulation, line emission, processes of formation and destruction, 
and mechanisms of population and depopulation of the molecu- 
lar energy levels as a function of the distance from the ionizing 
source of an ionized gaseous nebula. 

Although there is observational evidence that at least part of 
the H2 1-0 S(l) line emission may originate inside the ionized 
region of PNe, the published studies of the molecular hydrogen 
emission of PNe do not usually take into account the contribu- 
tion of the ionized region to this emission. One of the important 
conclusions of this work is that the H2 emission of the ionized 
work can contribute significantly to the total emission observed 
in PNe, particularly for PNe with high 7V Comparison between 
the calculated and observed H2 1-0 S(l)/Bry ratio shows that the 
emission of the ionized region can be responsible for a substan- 
tial fraction of the total H2 emission in such cases. Therefore it is 
important to use a code where the neutral and the ionized region 
are taken into account self-consistently. 

In the ionized region, the H2 IR emission lines are produced 
predominantly in the TZ, where the atomic lines [N II], [O I], and 
[S II] are significantly produced. The partially ionized and warm 
gas in the TZ favors the formation and survival of H2 molecules, 
as well as its IR line emission. The temperature of the central 
star is an important factor for the H2 density and IR line inten- 
sity, since hotter stars produce more high-energy photons than 
colder stars. These photons can penetrate deep into the nebulae, 
producing the TZ. The 1-0 S(l) line intensity increases strongly 
with the increase in in our models. This result agrees with the 
correlation between th e detection of H 2 in PNe and their Zanstra 
temperatures found by Phillip! (l2006h . Furthermore, this can ex- 
plain why the detection of this H2 line is more common in bipo- 
lar PNe (Gatley's ru le), given that the se objects typically have 
higher r*. Although lReav et all (Q988) suggest that the correla- 
tion between the intensity of the 1-0 S(l) and [OI] lines in PNe 
is due to the existence of clumps, such correlation can also be 
naturally explained by both lines being produced in the TZ. 

The most intense H2 lines are emitted in the range 1 to 29 
jt/m, in the bands 0-0, 1-0, 1-1, 2-0, 2-1, 3-1, 4-1, 4-2, and 5-3. 
Several of these intense lines have already been detected in PNe. 
The 1-0 S(l) line is one of the more intense lines. The fraction 
of the 1-0 S(l) line to the total H2 IR line emission is around 1% 
(within a factor of less than 10) in all our models. 

Both collisions and UV pumping play important roles in the 
excitation of H2 infrared lines in the ionized region. The effect 
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of the excitation by UV pumping is important for levels with 
v > 3. The relative importance of collisions over UV pumping 
increases with the increase in and n H - Grain surface reac- 
tion and associative detachment may only be significant for very 
high levels (v > 5 and J > 6), particularly in the cases of denser 
nebula or colder central stars. Grain surface reaction may be im- 
portant for lower levels in the presence of large amounts of dust. 

The effect of H2 on the thermal equilibrium is insignificant in 
most models, except in the cases of a very high central star tem- 
perature or dust density, where the molecule cools the gas in the 
outer zone of the TZ, primarily through collisional excitation. 
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Fig. 7. Ratios of some H2 lines to Bry for the ionized region as a function of . The models are the same as for Fig. [2] The observed 
ratios are represented by dots. Boxes with down arrows are upper limits. The references for the observations are given in the text. 
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